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Abstract 

For a discrete two-state quantum walk on the half-line, we formu- 
late and prove a weak limit theorem describing the terminal behavior 
of its transition probabilities. In this context, localization is possible 
even for a walk predicated on the assumption of homogeneity. For 
the Hadamard walk on the half line, the weak limit is shown to be 
independent of the initial coin state and to exhibit no localization. 

1 Introduction 

The notion of quantum walks (QW) originated in 1993 with the publication 
of p] . Since then, inspired by the ideas in [21 [3j H] , and propelled, in no small 
measure, by the hypothetical prospect of application to the development of 
super-efficient quantum algorithms, this field of research has flourished im- 
mensely. For a lively and informative elaboration of the history of quantum 
walks and their connection to quantum computing, the reader is referred to 

GaisiEiisiEinni. 

Among numerous other valuable results, Konno et al. [TT| lT2| FP3] have 
formulated and proved a weak limit theorem for the transition probabilities 
of discrete unitary quantum walks on the full line (dynamically restricted 
to integer nodes). In this work, we formulate and prove an analogous result 
for quantum walks on the half-line (dynamically restricted to non-negative 
integer nodes). Specifically, in this paper, we adopt the model of unitary 
discrete quantum walks on the half-line as proposed by Cantero et al. [H]. 

For a QW of this kind, much effort has been dedicated to studying the 
asymptotic behavior of its site-specific return probabilities and the mass 
point aspect of these distributions, notably by Cantero et al. using the 
CGMV method [TQl US]- Meanwhile, for a QW of this kind, the overall weak 
limit law for the probability distribution has remained relatively unexplored. 
Our treatment addresses this omission and emulates the approach employed 
by Chisaki et al. in |16j . 

This paper unfolds as follows. In the second section, we introduce the 
elements essential to understanding our adopted model of quantum walks 
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on the half line and we define and develop the generating functions of the 
probability amplitudes associated therewith. In third section, we formulate 
and prove a weak limit formula for the class of QW's under consideration. 
Much of the derivations and justifications of technical details are relegated 
to the subsequent sections. 



2 QWs on the half line 

In this paper, we adopt the model of unitary quantum walks on the half-line 
as defined by Cantero et al. [H]. In this model, the unitary shift operator 
S is given by 



s = j2\x + i)(x\®\ txt i + J2 \ x - ® i ^ 



®it>u 



2 = 1 



According to this formula, each temporal iteration of S acts as follows: 
when located at any site 0, 1, 2, 3, an upward spin migrates one step to the 
right. When located at any strictly positive position 1,2,3, a downward 
spin migrates one step to the left. When located at the special site 0, a 
downward spin reverses orientation and remains at site 0, so that, in the 
next iteration, it migrates to site 1. 

In this model, the evolution of quantum states is governed by two unitary 
so-called "coins" , represented by the matrices: 



U 



a 

c 



b 
d 



a b 
c d 

in terms of which, the one-step transition rules can be framed by a master 
operator it, as follows. When i ^ 0, we have 



U\xl) = a\x- 1 1) +c\x + 1 1), il|zt> = b\x-li) + d\x + l T), (3) 
while, in the special case when x = 0, we have 



il|0|) = a |0t)+c |lt), H|0T> = 6 |0 f> + d |l t>- (4) 

For the sake of clarity, it is appropriate to distinguish between two kinds 
of quantum walks on the half line. If the quantum walk is governed every- 
where by a single constant coin U (in particular Uq = U in Eq.Q) let it be 
referred to as a homogeneous quantum walk. Otherwise, let the quantum 
walk be called inhomogeneous. For instance, the Hadamard quantum walk 
on the half line is the homogeneous quantum walk governed everywhere by 
the constant coin U = H, where H is Hadamard matrix: 



As is customary in the literature, the probability amplitude of finding 
the walker at position x at time t is denoted by ip(x, t) = (ipi(x, t),ip^(x,t)) T , 
while its generating function is denoted by ^(x,z) = YltLo ^) z< - 

Referring to Eq. let the matrices 5*0, Qo, P and Q be defined as 
follows: 



So 



" 


, Qo — 


" 





, P = 


a b 


, Q = 


"00" 


_a Q b _ 




d _ 





c d 



(6) 



Then the one-step transition rules given by Eqs. @ and Q can be 
reformulated as follows in terms of the probability amplitudes: 



+ = Pt/>(l,t) + S i{>(0,t) 

V(l,t + 1) = P^(2,t) + Qo^(0,t) 

V»(a:,* + 1) = P^(x + l,t) + forx > 2 (7) 

Enlightened by the techniques employed by Konno and Oka et. al. in 
[TT| [T7] . it is possible, in the present context, to derive an explicit closed 
formula for the generating function *&(x, z). 

Lemma 1 Suppose the quantum walk (as defined above) is launched 
from site with initial coin state (a, (3) T (i.e., ^(0, 0) = (a, (3) T and if)(x, 0) = 
for x > 1). Then the generating function ^(x, z) is given by 

T , n \ (b + adB r (0,z))(aA z + aco + (3d )z 2 

^i{0,z) = a + 



* t (0,«) = p + 



l-b z- bc z 2 - bA z 3 - adA z 3 B r (0, z) - ac dz 2 B r (0, z) 

aa z + /3b z + (3bA z 3 + (3adA z 3 B r (0, z) 
l-b z- bc z 2 - bA z 3 - adA z 3 B r (0, z) - ac dz 2 B r (Q, z) 



T , x (d\ + /a) x 1 dzB r (0, z)(aA z + ac + /3d ) 

U ' J l-b z- bc z 2 - bA z 3 - adA z 3 B r {0, z) - ac dz 2 B r {0, z) ~ ' 

«, (r z) = (dX+fay-'ziaAoz + ac + (3d ) 

n ' J 1 -6 ^-^ 2 -6A z 3 -adA ^ 3 ^(0, aco^E^O, z) ~ ' 

where 

l + A Z 2 - v /A¥ + 2A(l-2|a|V + I Rrfn , A + - 

A+ " 2^ (8) 

A = a c? — b c , and A = ad — be (9) 

The details involved in deriving the above formulas for ty(x,z) are de- 
ferred to section 4. 

Let X t denote the position of the walker at time t. According to stan- 
dard quantum mechanical conventions, the probability distribution of X t is 
identified with \ip(x, t)\ 2 . In the sequel, equipped with the explicit expres- 
sions for *&(x, z), as given by Lemma 1, our main goal is to evaluate, for a 
QW on the half-line, the terminal behavior, in the weak limit sense, of the 
scaled quantity X t /t as t — > oo. 



3 Weak limit for QW's on the half line 



To enable feasibility within reasonable effort, and without excessive loss of 
generality, let us proceed on the assumption that the determinants of the 
unitary coin operators in Eq. ^ are equal. That is: Ao = A. Following 
the approach in [16J, the following theorem is obtained: 

Theorem 1. Suppose the QW on the half line (as defined above) is 
launched from site with the initial coin state a\ \) + (3\ f), then 

X t /t Y, (10) 

where Y is the random variable with probability density 

/H = *) + %)J%^= (ii) 
n 1 - y >V\ a \ - v 

and where the symbol =>• denotes convergence in distribution. The symbols 
in the expression for f(y) are defined as follows: 

(h & - ft-sX^i + h 2 - h 3 ) + h 7 (h 4 + h 5 ) 

h{y) = - 



(h + h 2 - h 3 ) 2 - (h 4 + h 5 ; 



2 



{he + hs){-hi + h 2 + h 3 ) + h 7 (h 4 - h 5 ) 

(-h 1 + h 2 + h 3 ) 2 -(h 4 -h 5 y ' 1 J 



where 



h = 2|c| 2 3(cAi - c Al), h 2 = (1 + |c | 2 - 2$t(cc ))\c\ v / l-y 2 , 
h 3 = 3(cA^ + cc 2 A-i)(l - y 2 ), h 4 = 2|c|3?(cA^ - cqA*) ^\a\ 2 - y 2 , 
h 5 = 2%cc Q )^l-y 2 ^\a\ 2 -y 2 , h e = (\a\ 2 + \ac + f3d \ 2 )^l-y 2 

h 7 = 23fJ(ctA5ac + f3d )^/\a\ 2 - y 2 , h 8 = 2\c\Q{aA^ac + f3d ), (13) 



and where 



\a\ 



h(y) |c|V/(0 ^ » (14) 

^(i -y 2 )V\ a \ -y 2 



The proof of this theorem can be found in Section 5. 

In [H1[I2], Konno offers weak limit theorems for homogeneous quantum 
walks on the full line. More recently, these results have been extended to 
topological settings similar [13], but not quite equivalent, to the half-line 
setting of Theorem 1. In our setting, when restricted to the special case 
of a homogeneous quantum walk on the half line, Theorem 1 implies the 
following corollary: 

Corollary 1. Let the homogeneous QW be launched from site with 
the initial coin state a\ l) + /3| t}- Then 



X t /t => Y, 



(15) 



where Y is the random variable with probability density 



r/ . 2|c| 3 [|a| 2 + \ac + (3d\ 2 - 2%(cA^)%(aA^ac + (3d)]y 2 T 

f{y) = (o(-A±\\2a /, | 2 2 J (o,H)U/) 

7r|a| 2 [|c| 2 — (b(cA2j) 2 (l — y 2 )\(l — 2r)y |a| — y 

In all cases studied by Konno in [ITJ [12], involving a two-state homoge- 
neous QW on the full line, the phenomenon of localization never occurs as 
a feature of the weak limit. However, for a typical two-state homogeneous 
QW on the half-line, the limiting behavior is quite different. Frequently, 
but not always, localization is to be expected. However, this discovery is 
not without precedent. Notably, Cantero et al. [10j[T5] have devised special 
techniques known as the CGMV method for treating questions about mass 
point and the return probabilities of QWs on the half line. Also they (see 
page 1159 in [11)]) have offered an example of a homogeneous QW on the 
half line for which localization occurs. Here we offer a different example. 

Example 1. Suppose the QW, everywhere governed by the coin oper- 
ator 



1 e^ /4 

e -i7r/4 _j 



(16) 



is launched from site with initial coin state ^| i) + T)- Then, by 
Corollary 1, we have 

f(y) = pS(y) + — , (17) 

n 1 -y Jv 1 - 2 v 

where p = (^ft^l . 

The coefficient p (see Eq. (14)) acts as the weight of the Dirac func- 
tion S(y) in the density function f(y) and thereby governs the propensity 
for localization near the origin. In fact it is not difficult to see that p may 
be computed as the sum of the limiting probabilities at all nodes, namely 
p = J2°?=o li m t->oo l^i^jt)] 2 . Accordingly, for the purposes of numerical sim- 
ulation, it suffices to consider only probability values not very far removed 
from the initial position and time, say x < 4 and t = 400, as in Table 1: 



T able 1: Probability Distribution of X t for x < 4 and t = 4 00 
X 1 2 3 I 

\ip(x, t) | 2 .0492 .0132 .0035 .0010 .0002 

In terms of the tabulated values, we get p ~ Y2 x =o lV'( a; '^)| 2 — 0.0671 ~ 
(V3-V2)(3-V3) ^ rjj^ ]j m j^j n g behavior of this specific QW is displayed in Fig- 
ure 1. Note the localization spike near the origin. 

However, as shown by our next example, localization is not always present 
for a homogeneous unitary QW on the half line. Let the Hadamard matrix 



_ ^^..aa^ 



Fi gure 1: Probability distribution at time t = 400 of a homogeneous QW on the 
half-line with unitary coin operator U as defined by Eq. (16 1 and initial coin state 

f U)+*^|t>- 



H (see Eq.([5])) serve as the constant coin operator for the homogeneous QW. 
By Corollary 1, we deduce: 

Corollary 2. Suppose the Hadamard QW is launched from site with 
the initial coin state a\ l) + j3\ f), then 



X t /t => Y 

where Y is a random variable with probability density 



;is) 



21 



(0,V2/2) 



(v) 



f(y) = 

7r(l-rVi_ 2 ^ 
For this QW, the probability distribution is depicted in Figure 2. 




Figure 2: Probability distribution at time t = 400 of a homogeneous QW on the half-line 
with Hadamard coin operator H as defined by Eq.M and initial coin state ^| ^ | t)- 



Notably, according to Corollary 2, at least when the QW on the half-line 
is governed by the Hadamard coin operator, the weak limit is seen to be 
independent of the initial coin state. 

Last, but not least, let us consider an example of an inhomogeneous 
quantum walk on the half line. 

Example 2. Let the QW on the half line be launched with initial coin 
state | f) and coin operators U and U defined as follows: 



Un 



V2 
2 



1 

e -j7r/4 



3 i7r/4 



U = H 



V2 
2 



1 1 
1 -1 



(19) 



By Theorem 1, the density for this QW is given by 



f(y) = Po%) 



hy 2 + k 2 y 4 + hy e 



'(0.V2/2) 



(y) 



(20) 



m 1 + m 2 y 2 + m 3 y 4 + m 4 y 6 ) vta/1 - 2y 2 ' 

where the numerical coefficients are 

p = 0.677887, h = 54 - 38^, k 2 = 46^ - 62, k 3 = 16V2 - 16 
mi = 17 - 12V2, m 2 = 64^ - 90, m 3 = 121 - 84^, m 4 = (4 - Ay/2) 2 . 

The relatively large value of p ~ 0.678 indicates a strong tendency for 
localization near the origin. In Table 2 and Figure 3 we display numerical 
data for this QW at time t = 400. 

Table 2: Probability Distribution of X t for x < 7 and t = 400 



x 


1 


2 


3 


4 


5 


6 


7 


\ip(x,t)\ 2 .4428 


.1648 


.0299 


.0275 


.0046 


.0054 


.0010 


.0012 



Using the tabulated data, we verify that po ~ XL=o l^/H 2 ^)! 2 — 0.6771. 




Figure 3: Probability distribution at time t = 400 of a non- homogeneous QW on the 



half-line with coin operators Uq and U as defined by Eq. ( 19 ) and initial coin state | t) 



Arguably, in the literature, the publication most closely related to our 
present treatment of weak limits of QWs on the half-line is the investigation 
by Chisaki et al. [18J of a QW on a graph consisting of joined half lines. 
In that work, a similar weak limit is obtained. Other relevant publications 

include [ini [201 1221 1231 [2H 123 

Recently, Konno et al. [13] have proposed a universal formula charac- 
terizing weak limit measures of quantum walks on the full line. Accord- 
ing to this formula, any QW on the line may be characterized as follows: 
there exists a rational polynomial w(y) and real numbers C G [0, 1) and 
r G (0, 1) such that the weak limit p(dy) = C5 (y) + w(y)fK(y', r)dy, where 



M^ r ) = n A -rMv) with C = 1- JZo w (y)fK(y;r)dy. We 

y )y r y 

remark that our QW model on the half line also conforms to this character- 
ization, with some adjustments. In particular, the indicator function I(- r ,r) 
over the symmetric interval (— r, r)would have to be replaced by I(o,r), where 
r = \a\. 



4 Proof of Lemma 1 



To obtain the desired generating functions, we begin by introducing some 
convenient notations. Let S(x, t) denote the transition amplitude of the 
given QW, defined as the sum of the contributions from all paths starting 
at the site and ending at the site x after t steps. For the specific model 
of a QW with a breakdown at the site 0, as in 

PH, let E b (x,t) denote 

3 6 (0 — > x,t) as in [17J. To denote the generating functions B r (0 —> x,z), 
B q (0 — > x,z) and B r (0 }" X) Zj 5 clS defined in [T7j, the abbreviated symbols 
£T(0, z), B q (x,z) and B r (x,z) respectively are introduced. 

In addition to the special matrices 5*0, Qo, P and Q defined in Eq. ([6]), 

let 



a b 




Rn 



c d 




s 





a b 



R 



c d 




(21) 



By the PQRS method [IT] , the transition amplitudes can be expressed 
as unique linear combinations of the matrices P , Qo, Rq and S (or P,Q,R 
and S): 



S(x, t) = u po (x, t)P + u qo (x, t)Qo + u ro (x, t)R + u so (x, t)S . 



(22) 



The expressions in Eq. (|22j) can be represented as coefficients of the 
following four generating functions in z: 



U po (x, z) = uP ° (x> uqo fa z) = uQ ° ( x > ( 23 ) 

t=0 t=Q 

oo oo 

U ro (x,z) = ^u ro (x,t)z 4 , U S0 (x,z) = ^u S0 (x 1 t)z t . (24) 
t=o t=o 

Now by Eq. ([7]), the transition amplitudes at x — 1 are given by: 

2(1,1) = Q 

"(1)2) = Qo$o — doSo 

E(l,t) = ^S 6 (0,t-l-/)g 2(0,/) 

i 

+ 3(0, t - 1)Q + Q 2(0, t - 1) for t > 3 (25) 

in terms of which, we obtain the following formulas for the generating 
functions at x — 1. 



U Po (l,z) = [c o U Po (0,z) + doU So (0,z)}dzB r (0,z), 

U qo (l,z) = z + d zU qo (0,z) + c zU ro (0,z), 

U r °(l,z) = [l + d o U qQ {0,z) + c o U ro {0,z)}dzB r {0,z), 

U S0 (l,z) = c o zU Po (0,z) + d zU s °(0,z). (26) 

Again by Eq. (]7]), we have: 



P5(l,t) + S E(p,t) = E(p,t + l) 



(27) 



Also, by Eqs. (|27j) and (26) we obtain the following formulas for the 
generating functions at x = 0: 



UP0 ^ Q s _ bd z 3 + add z 3 B r (0,z 



U go (0,z) 
U ro (0,z) 



l-b z- bc z 2 - bA z 3 - adA z 3 B r (0, z) - ac dz 2 B r (0, z) 

a bz 3 + aa dz 3 B r (0, z) 
l-b z- bc z 2 - bA z 3 - adA Q z 3 B r (0, z) - ac dz 2 B r (0, z) 

(1 -b z)[bz 2 + adz 2 B r (0,z)} 

l-b z- bc z 2 - bA z 3 - adA z 3 B r (0, z) - ac dz 2 B r (Q, z) 



rrsorn z) = z-bc z 3 -ac dz 3 B r (0,z) 

1 ' ; l-b z- bc z 2 - bA z 3 - adA z 3 B r (0, z) - ac dz 2 B r (0, z) ' 

When x > 2, Eq. @ yields: 



E(x,x) = E b (x — 1, x — l)Qo 

E(x,t) = £ 5 fe (x-l,i-l-/)Q H(0,0 

0<l<t-x 

+ E b (x -l,t- 1)Q for t > x + 1, (29) 

and subsequently: 



U P0 {x,z) = [coU Po (0,z) + d o U So (0,z)}dzB r (x-l,z), 

U qo (x lZ ) = [l + d U qo (0,z) + c U ro (0,z)}dzB q (x~l } z) } 

U ro (x lZ ) = [l + d U qo (0,z) + c U ro (0,z)}dzB r (x-l } z) } 

U s °(x,z) = [c o U Po (0,z) + d U s °(0,z)]dzB q (x - l,z). (30) 



Finally, we have: 



V±(0,z) = (1,0)^(0,0)+ ^(1,0)5(1,^(0,0); 

t=i 

= a + (ac + f3d )U r °(0,z) + (aa + pb o )U Po (0, z) 

oo 

%(x,z) = ^(l,0)H(M)V(0,oy 

t=\ 

= (ac + /3d )U r °(x, z) + (aa + /3b )U P0 (x, z) 

tf t (0,z) = (0,1)#,0) + ^(0, l)E(x,t)^(0,0)z t 

t=i 

= (3 + (ac + (3d )U qo (0, z) + (aa + f3b )U s °(0, z) 

oo 

^{x,z) = ^2(0,l)E(x,t)iP(0,0)z t 
t=\ 

= (ac Q + pd )U qo (x,z) + (aa Q + pb )U S0 (x,z) (31) 



Together, Eqs.fl26b, (28), l30h and (31) suffice to complete the proof. 



5 Proof of Theorem 1 

Consider the Fourier transform of the generating function $?(x,z), which 
is given by *&(k, z) = Y^Lo ^{ x i z)e lxk . Collecting like terms in z f , we get 
$(k,z) = Et=o^(M)^ where *(M) = Y.T=o^i. x ^) eixk - % Lemma 1, 
one has a closed formula for ^f(k, z), by which it can be shown that ty(k, z) 
has two types of unit poles: the first type of poles are the zeros (assumed to 
be {uj :j E J}) of l-b z-bc z 2 -bA z 3 - adA z 3 B r (0, z) -ac dz 2 B r (0, z), 
which contribute to the part of a mass point represented by the distribution 
pS(y) in the formula for f(y); while the second type of pole, as the zero of 
1 — ^A + e* fc , determines the absolutely continuous part of the density f(y). 
Let us introduce aA~^ = \a\e iri with < rj < 2ir, then this single pole 
is e t9l( - k ^ when k E (r) — 7r, 77) , or e %e ^ k ' when k E (j] — 2n,r] — it). Here 
9i(k) = arccos(|a| cos(?7 — k)) and 9 2 (k) = 2n — arccos(|a| cos(r? — k)). 
By Cauchy's integral theorem, 



„ n 1 / y(k,z) , . . 

*(M) = — f -^dz (32) 

Z7 ™ </-j>:|.2|=rfor0<r<l} z 

Res(^,e^) - E^Res^,^-), k E (rj — tt,^ 



-Res(™, e^) - V. eJ R es (™ ^), k E (rj — 2n, rj 



7T 



By the structure of the function \l/(/c, t), one can calculate the character- 



istic function of the scaled quantity X t /t so that 

dk 

27T 



E(e^) = / (9(k,t)Mk + Z/t,t)) : 

Jn-2n 



r)—2n 

7,-TT 27T 

/ e^ t+1 ^ k ^ k+ ^Res(^(k, z), e^( fe ))Res(#(A; + f/t, z), e w ^ k +^) — 

T V^R M k > Z ) , *{k + t/t,z) dk 

/ 2^ Res (^T^' w i) 2^ Res ( ^1 ' w i) 2^ 



+ e l(m)9l(fc) Res($(fc,2),e^( fc ))^Res( 

•/ 77— 7T „v- 7 



V(k + £/t,z) ^dk 

CO,; 



t +l ' 3>2ir 

>„-r. j£j 



+ [" e -(m)^A) Reg( ^ + ^ 

+ [ V ~* e^ +1 >^Res($(fc, z), e^W) V Res( ^ ffi z) , 



+ /" ' e- J (* +1 ) e2 ( fe+ ^Res($(A; + £/t, z), e *> (**/*)) V R es (*&^, (34) 



As t — > oo, each of the last four terms in Eq. (34) converges to by 
Riemann-Lebesgue lemma, and the third term converges to lim^oo J2T=o \^{ x i I 2 ' 
which is p in the formula for the density function f(y). Therefore we get 

, dk 



\imE(e^ Xt/t ) = f e-^ (fc) ||Res($(fc,2),e 
+ T V^ 2(fc) ||Res(§(£;,z),e' 

Jri-2-K 



ifl(fc)^||2' 

2n 
dk 



~ l ^ k) \\Res(^(k,z),e ie ^ k) )\\ 2 

J7-27T 27T 

+ P (35) 

Hsinfa-fc) j ol (U\ — |a|ginfa-fc) 



where ^(Jfc) = , MM an d ^(fc) 



Let y = | |a|sm(rj fcj ^ ^fter an extensive calculation, Eq. (|35|) re- 

•v/l — |a| 2 cos 2 (77— k) 



^Jl—\a\ 2 cos 2 (tj— k) -y/l— |a| 2 cos 2 (77— fc) 

|a| sin(' 
-|a| 2 o 

duces to the identity: 

lim E(e^) = p + r e^h(y)— ^ dy (36) 

5 "+ 00 Jo n(l - y 2 )J\a\ - y 2 



with h(y) is given by Eqs. (12) and (13), whereby the proof is complete. 
Acknowledgments 

Undergraduate research assistant and mathematics major Forrest Ingram- 
Johnson is to be commended for performing the computer simulations for 
the graphs in Fig. 1-3. Liu was supported, in part, by NSF Grant CCF- 
1005564. 



References 



[I] Y. Aharonov, L. Davidovich, and N. Zagury, Phys. Rev. A 48, 1687 
(1993). 

[2] A. Ambainis, E. Bach, A. Nayak, A. Vishwanath and J. Watrous, One- 
dimensional quantum walks, in Proceedings of the thirty-third annual 
ACM symposium on Theory of computing, STOC '01 (ACM, New York, 
NY, USA, 2001) pp. 37-49. 

[3] D. Aharanov, A. Ambainis, J. Kempe and U. Vazirani, Quantum walks 
on graphs, in Proceedings of the 33rd Annual ACM Symposium on The- 
ory of Computing, (ACM, New York, 2001), pp.50-59. 

[4] E. Farhi, S. Gutmann, Quantum computation and decision trees, Phys. 
Rev. A 58, 915928 (1998). 

[5] J. Kempe (2003), Quantum random walks - an introductory overview, 
Contemp. Phys. 44, 307. 

[6] V. Kendon, Decoherence in quantum walks - a review, Struct, in Comp. 
Sci 17(6) pp 1169-1220 (2006). 

[7] N. Konno, in Quantum Potential Theory, Lecture Notes in Mathemat- 
ics, edited by U. Franz and M. Schurmann (Springer- Verlag, Heidelberg, 
2008), pp.309-452. 

[8] S.E. Venegas-Andraca, Quantum Walks for Computer Scientists, Morgan 
and Claypool Publishers (Synthesis Lectures on Quantum Computing), 
2008. 

[9] S.E. Venegas-Andraca, Quantum walks: a comprehensive review, Quan- 
tum Information Processing vol. 11(5), pp. 1015-1106 (2012). 

[10] M. J. Cantero, F. A. Griinbaum, L. Moral and L. Velazquez (2012), 
The CGMV method for quantum wafe, Quantum Inf Process (2012) 
11:11491192. 

[II] N. Konno, Quantum random walks in one dimension, Quantum Infor- 
mation Processing, 1:, pp. 345-354 (2002). 

[12] N. Konno, A new type of limit theorems for the one- dimensional quan- 
tum random walk, Journal of the Mathematical Society of Japan, 57: 
1179-1195 (2005). 

[13] N. Konno, T. Luczak, and E. Segawa, Limit measures of inhomoge- 
neous discrete-time quantum walks in one dimension,Q,uantum Inf Pro- 
cess (2013) 12:3353. 



[14] M. J. Cantero, F. A. Griinbaum, L. Moral and L. Velazquez (2010), 
Matrix-valued Szego polynomials and quantum random walks, Communi- 
cations on Pure and Applied Mathematics, 63, 464-507. 

[15] M. J. Cantero, F. A. Griinbaum, L. Moral and L. Velazquez (2012), 
One- dimensional quantum walks with one defect. Rev. Math. Phys. 24, 
1250002. 

[16] K. Chisaki, M. Hamada, N. Konno, E. Segawa (2009), Limit theorems 
for discrete-time quantum walks on trees, Interdiscip. Inform. Sci. 15, 
423-429. 

[17] Oka, T., Konno, N., Arita, R., and Aoki, H., Breakdown of an electric- 
field driven system: a mapping to a quantum walk, Physical Review 
Letter, 94: 100602 (2005). 

[18] K. Chisaki, N. Konno, and E. Segawa, Limit Theorems for the Discrete- 
Time Quantum Walk on a Graph with Joined Half Lines, Quantum In- 
formation and Computation 12, (2012) 0314-0333. 

[19] G. Grimmett, S. Janson and P. F. Scudo (2004), Weak limits for quan- 
tum random walks, Phys. Rev. E 69, 026119. 

[20] N.Inui, N. Konno and E. Segawa, One- dimensional three-state quantum 
walk, Phys. Rev. E 72 (2005) 056112. 

[21] M. Katori, S. Fujino and N. Konno, Quantum walks and orbital states 
of a Weyl particle. Phys. Rev. A 72 (2005) 012316. 

[22] T. Miyazaki, M. Katori and N. Konno, Wigner formula of rotation 
matrices and quantum walks. Phys. Rev. A 76 (2007) 012332. 

[23] E. Segawa, N. Konno, Limit theorems for quantum walks driven by 
many coins. Int. J. Quantum Inf. 6 1231 (2008). 

[24] C. Liu, N. Petulante, One- dimensional quantum random walk with two 
entangled coins. Phy. Rev. A 79 032312 (2009). 

[25] N. Konno, T. Machida, Limit theorems for quantum walks with memory. 
Quant. Inf. Comput. 10 1004 (2010). 

[26] C. Liu, Asymptotic distributions of quantum walks on the line with two 
entangled coins, Quantum Inf Process (2012) 11:1193-1205. 



